Josephson coupling through ferromagnetic heterojunctions with non-collinear 
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We study the Josephson effect in clean heterojunctions that consist of superconductors connected 
through two metallic ferromagnets with insulating interfaces. We solve the scattering problem 
based on the Bogoliubov-de Gennes equation for any relative orientation of in-plane magnetizations, 
arbitrary transparency of interfaces, and mismatch of Fermi wave vectors. Both spin singlet and 
triplet superconducting correlations are taken into account, and the Josephson current is calculated 
as a function of the ferromagnetic layers thicknesses and of the angle a between their magnetizations. 
We find that the critical Josephson current I c is a monotonic function of a when the junction is 
far enough from — n transitions. This holds when ferromagnets are relatively weak. For stronger 
ferromagnets, variation of a induces switching between and n states and I c {ot) is non-monotonic 
function, displaying characteristic dips at the transitions. However, the non-monotonicity is the 
effect of a weaker influence of the exchange potential in the case of non-parallel magnetizations. No 
substantial impact of spin-triplet superconducting correlations on the Josephson current has been 
found in the clean limit. Experimental control of the critical current and — n transitions by varying 
the angle between magnetizations is suggested. 

PACS numbers: 74.50.+r, 74.45.+C 
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I. INTRODUCTION 

The interplay between ferromagnetism and supercon- 
ductivity in superconductor (S) - ferromagnet (F) hy- 
brid structures attracts considerable interest for some 
time alreadyi^* 2 * 3 ^ Variety of interesting theoretical pre- 
dictions, such as the existence of 7r-state superconductiv- 
ity in SF multilayers; 5 ^ 7 ^ 9 ! 10 ! 11 ! 12 ! 13 has been confirmed 
experimentally 14 i 15 i 16 i 17 i 18 i 19 i 20 In SFS Josephson junc- 
tions, the Zeeman effect induces both a strong decay and 
oscillations of superconducting correlations inside the 
ferromagnet The corresponding oscillations of the 
Josephson critical current with thickness of the ferromag- 
netic layer or the exchange energy have been calculated 
in both the clean and the dirty limit in the framework 
of quasiclassical theory of superconductivity! 7 ^^ Experi- 
mental evidence for such oscillating behavior in weak fer- 
romagnetic alloys came only recentlyii&iiiSSi Similar ex- 
periments with strongly spin-polarized ferromagnets 2 ^ 2 ^ 
have been performed in spintronic setup^ 5 - However, 
for strong ferromagnets and finite transparency of inter- 
faces the quasiclassical treatment is no longer valid and 
Bogoliubov-de Gennes or Gor'kov equations have to be 
solved^SJ. 

Heterostructures with superconductors coupled 
through inhomogeneous ferromagnets have been also 
extensively studiedi 28 i 29 i 30 i 31 i 32 i 33 i 34 i 35 The simplest 
example of such a structure contains a two-domain ferro- 
magnet with collinear magnetizations, either parallel (P) 
or antiparallel (AP). In this case, besides spin singlet su- 
perconducting correlations, only triplet correlations with 
zero total spin projection exist i^SiSLS These correlations 
penetrate into the ferromagnet over a short length scale 



determined by the exchange energy. For non-collinear 
magnetizations, triplet correlations with nonzero spin 
projection are present as well; they are not suppressed 
by the exchange energy, and consequently they are 
long-rangedi^ It is predicted that triplet components 
should have a dramatic impact on the Josephson effect, 
displayed through a non-monotonic dependence of the 
critical current on angle between magnetizations. 29 ' 30,31 
The triplet correlations were proposed as a possible 
explanation for recent observations of a Josephson 
current through half-metallic barriers^i 4 ^ In diffusive 
Josephson junctions, the length scales associated with 
short- and long-range correlations are, respectively, 
£/ = ^hDf/ho and £ s = ^JhjJJjk^T , where Df is 
the diffusion constant in the ferromagnet, and thermal 
energy k^T is typically much smaller than the exchange 
energy ho. It is desirable to know if such interesting 
effects are also relevant in ballistic SF heterostructures, 
where the ferromagnet coherence length £f = hVp/ho 

is the only characteristic length, v F being the Fermi 
velocity. 

In this paper we investigate the Josephson effect in a 
clean SI1F1I2F2I3S heterojunction where the ferromag- 
netic interlayer consists of two mono-domain layers hav- 
ing a relative angle a between their in-plane magneti- 
zations. Using the Bogoliubov-de Gennes formalism, we 
calculate the Josephson current and demonstrate that its 
critical value I c is a monotonic function of angle a when 
the junction is far enough from — n transitions. How- 
ever, this is possible when ferromagnets are weak. For 
stronger ferromagnets, I c (a) is a non-monotonic func- 
tion of a with characteristic dips related to the onset 
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of — 7r transitions. This non-monotonicity is a simple 
consequence of a weaker influence of the exchange po- 
tential in the case of non-parallel magnetizations. No 
substantial impact of spin-triplet superconducting corre- 
lations on the Josephson current has been found in the 
clean limit. In fully transparent symmetric junctions the 
ferromagnetic influence is practically cancelled out for 
antiparallel magnetizations. This is not the case for fi- 
nite transparency of interfaces when coherent geometrical 
oscillations of the critical Josephson current are superim- 
posed on oscillations related to the transitions between 
and 7r states. In the limiting case of zero exchange en- 
ergy our solutions reduce to the results for SI1NT2N2I3S 
junctions with complex non-ferromagnetic normal-metal 
(N) interlayer, previously studied in the quasiclassical 
approach 4i 

The paper is organized as follows. In Sec. II we present 
the model and solutions of the scattering problem used 
for calculation of the Josephson current. In Sec. Ill we 
discuss numerical results for the current-phase relation, 
and dependence of critical current on angle between mag- 
netizations and thickness of ferromagnetic layers. We il- 
lustrate these results for various spin polarizations and 
transparencies of the interfaces. Concluding remarks are 
given in Sec. IV. 



II. MODEL AND SOLUTION 

We consider a clean SI1F1I2F2I3S heterojunction con- 
sisting of superconductors (S), two uniform mono-domain 
ferromagnetic layers (Fi and F2) with misorientation an- 
gle a between their magnetizations, and nonmagnetic in- 
terfacial potential barriers between metallic layers (Il- 
ls), see Fig. ^ Superconductors are described in the 
framework of standard BCS formalism, while for ferro- 
magnets we use the Stoner model with an exchange en- 
ergy shift 2h (r) between the spin subbands. 
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FIG. 1: Schematics of an SI1F1I2F2I3S heterojunction. The 
magnetization vectors lie in the y-z plane and form the oppo- 
site angles ±a/2 with respect to the z-axis. 



Electron-like and hole-like quasiparticles with energy 
E and spin projection a =|,i are described by wave- 
functions u a {y) and v a (r), where r is the spatial coordi- 
nate. Using the four-component wave function \&(r) = 



[tt|(r), U|(r), «|(r), U|(r)] T , we write the Bogoliubov-de 
Gennes equation as 



where 



H 



H*(r) = E*(r), 



H+(r) A(r) 
A*(r) -F-(r) 



(2.1) 



(2.2) 



while the 2x2 blocks are given by H±(r) — Ho(r)l — 
h(r) sina(r)o : 2 T h(r) cos a(r)as and A(r) = A(r)l. 
Here, <7j and 1 are the Pauli and unity matrix, respec- 
tively, and H (r) = -h 2 V 2 /2m + W(r) + U(r) - /z. The 
chemical potential is denoted by [i, W(r) = J2i WiS(x — 
Xi) is potential of the barriers at interfaces, and U(r) is 
the electrostatic potential. The x-axis is chosen to be 
perpendicular to the layers, whereas x\ = —d\, x^ = 0, 
and X3 = e?2 are coordinates of the interfaces. The 
electron effective mass m is assumed to be the same 
throughout the layers. The difference [i—U(y) is equal to 

(s) 

the Fermi energy of superconductors, E F , or the mean 
Fermi energy of ferromagnets, E ( F f) = (El+E F )/2. Mod- 



uli of the Fermi wave vectors, k 



(s) 



2mE ( F >/h 2 and 



„(/) 



2mE F f] /h 2 in S and F may be different in gen- 
eral. The in-plane, y-z, magnetizations of the neighbor- 
ing F layers are not collinear in general, and the magnetic 
domain structure is described by the angle a(r) with re- 
spect to the z-axes: a(r) = —a/2 for — di < x < in Fi, 
and a(r) = a/2 for < x < g?2 in F2. We assume equal 
magnitudes of the exchange interaction in ferromagnetic 
domains, h(r) — ho<d(x+di)0(d2—x), where Q(x) stands 
for the Heaviside step function. We also assume that the 
two superconductors are identical, and neglecting self- 
consistency, we take the pair potential A(r) in the form 



A(r) 



di) 



2 0{x-d 2 



(2.3) 



where A is the bulk superconducting gap and cf> is 
the macroscopic phase difference across the junction. 
The temperature dependence of A is given by A(T) = 

A(0) tanh (l.74y/T c /T-l) 



Note that self-consistency may be safely neglected 
when the proximity effect is weak between S and F layers. 
This includes the situations with large tunnel barriers at 
interfaces, and/or narrow F constriction, and/or large 
Fermi velocities mismatch. For a planar junction geome- 
try and good contacts between S and F layers, A will be 
suppressed at the vicinity of FS interfaces. 

The parallel component of the wave vector, k||, is con- 
served due to translational invariance of the junction in 
directions perpendicular to the x-axis. Consequently, the 
wave function can be written in the form 



\Er(r) = ip(x)e 



(2.4) 
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where ip(x) = [itf(x), v±(x), v-\ (x)] T satisfies the 

boundary conditions 



ip(x)\xt+o = ip{x)\ Xi -o 



dip(x) 



dx 



dtp(x) 



Xi+o dx 

2lM 



Xi-0 



1>( Xi ), (2.5) 
Zik^ipfa). (2.6) 



Here, Zi = 2mWi/h k F ' (i = 1, 2, 3) arc parameters that 
measure the strength of each insulating interface located 
at Xi — —di, 0, d%. 

The four independent solutions of the scattering prob- 
lem for Eq. H2.1JI correspond to the four types of quasi- 
particle injection processes: an electron-like or a hole-like 
quasiparticle injected from either the left or from the 
right superconducting electrode^ When an electron-like 
quasiparticle is injected from the left superconductor, the 
solutions of Eq. (|2.1|l are 



u\{x) 


= (e ik+x 4 
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U[{x) 


= (e ik+x 4 


- bie 
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x ve~ 
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(2.8) 




= (e ik+x 4 


- b-\e 


~ lk+x )v^ 


- a-\e ik ~ x u, 


(2.9) 


Vf(x) 


= (e ik+x 4 


- h e ~ 


~ lk+x )v^ 


- ai e ik ~ x u, 


(2.10) 



for the left superconductor (x < 0), 



u T (x) 


= c T e ifc+a; we^ /2 4- 


die- lk 




= c x e ifc+:c we^ /2 4- 


die~ lk 




= c T e ifc+a: i; 4- d;e~ 


~ lk ~ x u, 


w T (x) 


= c l e ik+x v + d^ 


~ lk ~ x u, 



' x ve^ 2 , 



J^>/2 



(2.11) 
(2.12) 
(2.13) 
(2.14) 

for the right superconductor (x > d%), and 

u T (a?) = Cie<"cos(a/2)4-C , 2 e"< a; cos(a/2) 

-iC 3 e lq t x sin (a/2) - ide - * 8 * * sin (a/2), 

(2.15) 

U; (x) = ~iC ie lq t x sin(a/2) - iC 2 e- iq t x sm(a/2) 

+C 3 e ig i x cos (a/2) 4- CaT^t* cos (a/2), 

(2.16) 

Wi (a;) = iC 5 e lq i x sin (a / 2) +iC 6 e-* q i x sin (a/ 2) 

+C 7 e lq i x cos (a/2) 4- C 8 e~ lq i x cos (a/2), 

(2.17) 

v r (x) = C 5 e lq T x cos(a/2)+C 6 e-* q i x cos(a/2) 

+iC 7 e zq i x sin (a/2) + iC^e^i' sin (a/2), 

(2.18) 

for the left ferromagnetic layer Fi (— d\ < x < 0). Solu- 
tions for the right ferromagnetic layer F 2 (0 < x < d 2 ) 
can be obtain by substitution a — > —a, with a new 



set of constants C' X ,...,C' % . Here, u = + fl/E)/2 
and w = — Q/E)/2 are the BCS amplitudes, and 
fi = \/E 2 — A 2 . Longitudinal (x-) components of the 
wave vectors are 



fc± = s j{2m/K 2 ){E F s) ± fi) - k 2 



in superconductors, and 



<?± = sJ(2m/h 2 )(E F f) + p a h ±E)-k 2 



(2.19) 



(2.20) 



in ferromagnetic layers. The sign 4- or — in the super- 
script corresponds to the sign of quasiparticle energy, 
whereas p a = 4-l(— 1) is related to the spin projection 

u = T(i)- 

Note that solutions for quasiparticlcs with opposite 
spin orientation in F layers are coupled due to misorienta- 
tion of magnetizations. Solutions are decoupled for a = 
and a = 7r. In these cases spin triplet superconducting 
correlations with nonzero spin projections also vanish. 
The case of transparent SFFS junctions with a — tt was 
treated previously by Blanter and Hekkingi^l For a = 0, 
Z<i = 0, and Z\ — Z 3 solutions reduce to analytic expres- 
sions previously obtained for SIFIS junctions 2 ^! with a 
net thickness of the ferromagnetic layer d — d\ 4- in- 
complete solution of the scattering problem requires 
determination of 24 unknown coefficients: 44-4 constants 
in superconducting electrodes [a a and b a in the left S, 
and c a and d a in the right S for two spin orientations) and 
8 4-8 constants in ferromagnetic electrodes (Ci, . . . ,Cg 
in Fi and C 1; . . . , C 8 in F2). When applied to the so- 
lutions given by Eqs. ^577 ^ — ^5718^ . the boundary condi- 
tions at three interfaces, Eq. I|2.5|l . provide the necessary 
24 equations. Analogously, one can find the solutions for 
other three channels of quasiparticle injection processes. 
However, the Andreev amplitudes in the first channel are 
sufficient to calculate the Josephson current^ 



III. JOSEPHSON CURRENT 

The stationary Josephson current can be expressed in 
terms of the Andreev reflection amplitudes, a a = a CT (</>), 
by using the temperature Green's function formalism^ 



tr.ki 



^2K {k " +K - 



(3.1) 



where fc+, k~, and a an (<p) are obtained from k~ , 
and a a {_4>) by the analytic continuation E — * iw n . The 
Matsubara frequencies are lu ti = ■nkBT{2n 4- 1), n — 
0, ±1, ±2, and fi„ = y/uj^ + A 2 . 

For simplicity, we illustrate our results for a single 
transverse channel, k|| = 0, symmetric junctions, d\ = 



4 



e?2 = d/2, and equal Fermi wave vectors, kp = k F = 
hp. Superconductors are characterized by fixed bulk 
value of the pair potential, A/Ep^ 1 = 10~ 3 . In order 
to adhere to realistic parameters all the examples are 
shown for kpd < 100 since it still seems unlikely that the 
ballistic limit can be experimentally achieved for thicker 
ferromagnets. 
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FIG. 4: (Color online) The current-phase relation I(<j>) for 
T/T c = 0.1, Ii /Ef = 0.1, Zi = Z 2 = Z 3 = 0, dk F = 30, 
and three different values of the misorientation angle: a = 
(solid curve), a = n/2 (dashed curve), and a — n (dash- 
dotted curve). The latter approximately coincides with I c {d) 
for the corresponding SNS junction (dotted curve). 



FIG. 2: (Color online) The current-phase relation I(cj>) for 
T/T c = 0.1, h /E F = 0.01, Zi = Zi = Z 3 = 0, dk F = 30, and 
three different values of the misorientation angle: a = (solid 
curve), a = n/2 (dashed curve), and a — n (dash-dotted 
curve). The latter coincides with I ((f)) for the corresponding 
SNS junction (h = 0). 
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FIG. 5: (Color online) The critical current I c as a function 
of d for T/T c = 0.1, h /E F = 0.1, Z x = Z 2 = Z 3 = 0, and 
three different values of the misorientation angle: a = (solid 
curve), a = n/2 (dashed curve), and a — tt (dash-dotted 
curve). Dotted curve represents I c (d) for the corresponding 
SNS junction (h = 0). 



FIG. 3: (Color online) The critical current I c as a function 
of d for T/T c = 0.1, h /E F = 0.01, Z x = Z 2 = Z 3 = 0, 
and three different values of the misorientation angle: a = 
(solid curve), a = n/2 (dashed curve), and a — n (dash- 
dotted curve). The latter approximately coincides with I c (d) 
for the corresponding SNS junction (dotted curve). Inset: I c 
as a function of a for dk F — 30. 



For fully transparent interfaces, Z\ = Z2 = Zj, = 0, 
and for a weak exchange field, ho/Ep = 0.01, the current- 
phase relations are shown in Fig. [21 for three different 



values of misorientation angle a — 0, 7r/2, tt. The figure 
illustrates a junction with thin ferromagnetic interlayer, 
k F d = 30, at low temperature T/T c = 0.1. As expected, 
the I(4>) curve for a — tt is the same as for the corre- 
sponding SNS junction (ho = 0) because the influence of 
opposite magnetizations in F\ and F% practically cancels 
out. The I(4>) relations are anharmonic due to full trans- 
parency and low temperature. The junction is clearly in 
the 0-state for all values of misorientation angle; for such 
a small ho there are no oscillations of the critical current 
related to the — tt transitions for kpd < 100, Fig. [3] 
We have shown that the critical current increases mono- 
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FIG. 6: (Color online) The critical current I c as a function of 
the misorientation angle a for T/T c = 0.1, ho/Ep = 0.1, and 
four sets of interface transparencies: Zi = Z 2 ~ Zg — (solid 
curve), Z\ — Z$ — 0, Z2 = 1 (dash-dotted curve), Z\ — Zz — 
1, Z2 — (dashed curve), and Z\ — Z 2 = Zz = 1 (dotted 
curve). Panel (a): dkp = 30. Panel (b): dkp = 60. The dips 
in I c (pt) separate alternating and 7r states, beginning with 
a state from the left. 



FIG. 7: (Color online) The critical current I c as a function 
of d for an SFIFS junction with h /E F = 0.1 at T/T c = 0.1, 
and three different values of the misorientation angle: a = 
(solid curve), a = n/2 (dashed curve), and a — tt (dash- 
dotted curve). Panel (a): Z\ = Z 3 = 0, Z 2 = 1. Panel (b): 
Z x = Z 3 = 0, Z 2 = 10. I c (d) for the corresponding SNINS 
junction (dotted curves) are shown for comparison. 



tonicaly with increasing a, from the SFS (a — 0) to the 
SNS (a — 7r) value. It is evident that spin triplet corre- 
lations, present for < a < n, here do not induce a non- 
monotonic variation of I c (a). For larger exchange en- 
ergy, ho/Ep = 0.1, the current-phase relations are shown 
in Fig. 0] for the same values of a, kpd, and T/T c as in 
Fig. El In Fig. for a = ir/2 the I(<f) curve exhibits co- 
existence of and 7r states with a dominant contribution 
of the second harmonic^ The critical Josephson current 
oscillates as a function of the ferromagnetic interlayer 
thickness due to the onset of — 7r transitions, except for 
a close to 7T. The period of oscillations increases with 
a. This is simply the effect of a weaker influence of the 
exchange potential in the case of non-parallel magneti- 
zations. For example, the period of oscillations in I c is 
approximately twice bigger for a = tt/2 than for a — 0. 
We have obtained the same result for parallel, but twice 
smaller, exchange energy (a = and 1i /Ef = 0.05). It 
can be shown (from the calculation of the junction free 
energy) that the dips in I c {d) curves correspond to — tt 
transitions. Similar oscillations of the critical Josephson 
current as a function of a for fixed k-pd (either 30 or 60) 
are shown in Fig. [SI as well as the influence of finite in- 
terfacial transparency. As already suggested for diffusive 
junctions^i^ it can be clearly seen that the transition 
between and tt states can be induced by varying the 
misorientation angle a. This feature is of particular im- 
portance for experimental applications, since fine tuning 



of — 7r transitions could be realized more easily by vary- 
ing a than by varying thickness or temperature. 

The influence of geometrical resonances is illustrated 
in Fig. □ for Z x = Z z = 0, and (a) Z 2 = 1 and (b) 
Z2 = 10. All other parameters are the same as in Fig. [5] 
Oscillations of the critical current due to geometrical res- 
onances are superimposed on the oscillations related to 
the transitions between and 7r states. This effect is 
clearly visible in Fig. \7\b) for low transparency of the 
interface between magnetic domains (Z2 = 10). We em- 
phasize that now, due to the presence of a barrier be- 
tween the F layers, the magnetic influence persists for 
a = tt, and I c (d) considerably differs from the results ob- 
tained for the corresponding SNINS junction. However, 
for identical domains, there are still no transitions to the 
7T state in the AP configuration. It can be seen that the 
amplitudes of geometrical oscillations of the supercurrent 
are significantly larger in SFIFS than in the correspond- 
ing SNINS junctions (dotted curves in Fig. [7J. Rapid 
oscillation of I c {d) can be also seen in the case of finite 
transparency of interfaces between ferromagnetic layers 
and superconductors, 2 ^ for any value of a. This is due 
to the resonant amplification of the Josephson current 
by quasi-bound states entering the superconducting gap 
as the thickness d of the interlayer is varied4^ The ef- 
fect of finite transparency of all the interfaces is similar. 
In a more realistic case of planar (multichannel) junc- 
tions these geometrical oscillations of the critical current 
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are damped, while positions of maxima and minima are 
slightly shifted^ 

Transitions between and tt states can be induced by 
changing the temperature of a junction with low trans- 
parency and strong ferromagnetic influence. 26-46,47 This 
is illustrated in Fig. |S| for Z\ = Z% = 0, and Z2 = 10, 
k-pd = 43, ho/Ep = 0.1 and for three values of a = 0, 
7r/16, 7r/2. It can be seen that temperature dependence 
is very sensitive to the value of a; for example, the tem- 
perature induced transition from to tt state when a ~ 
is absent for a = ir/2. 
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FIG. 8: (Color online) The critical current I c as a function 
of T/T c for dk F = 43, h /E F = 0.1, Zi = Z A = 0, Z 2 = 10 
and three different values of the misorientation angle: a = 
(solid curve), a — 7r/16 (dash-dotted curve), and a — n/2 
(dashed curve). For a ~ the characteristic nonmonotonic 
variation is related to the transition from state (low T) to 
7r state (higher T) with dip at the transition. For a — n/2 
the junction is in state at all temperatures. 



two mono-domain ferromagnetic layers with arbitrary 
transparency of the interfaces and any angle a between 
magnetizations, including two previously considered lim- 
iting cases a — and irm££^ The Josephson current is 
calculated numerically via Bogoliubov-de Gcnnes formal- 
ism. We have found that transitions between and tt 
states, resulting in characteristic dips in I c (a) curves, can 
be induced by varying the relative orientation of magne- 
tizations, like in diffusive j unctions However, this is 
simply the effect of a weaker influence of the exchange 
potential in the case of non-parallel magnetizations. No 
substantial impact of spin-triplet superconducting cor- 
relations on the Josephson current has been found in 
the clean limit. For weak ferromagnets, far from — tt 
transitions, the critical Josephson current monotonically 
depends on the angle between magnetizations. While 
in fully transparent junctions oscillatory dependence of 
the critical Josephson current on junction parameters is 
related only to the — tt transitions, for finite trans- 
parency of interfaces pronounced geometrical oscillations 
occur due to coherent contribution of quasiparticle trans- 
mission resonances (quasi-bound states) to the Andreev 
process. 

In conclusion, due to recent progress in nanofabrication 
techniques^ the SI1F1F2I2S junction may be realized in 
the clean regime in a setup where the angle a could be 
tuned by a weak external magnetic field. We thus sug- 
gest experimental investigation of the predicted control 
of the critical current and — tt transitions by varying 
the angle between magnetizations in Josephson junctions 
with magnetic bilayers. 



IV. CONCLUSION 

We have studied the Josephson effect in clean 
superconductor-ferromagnet heterojunctions containing 
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